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Abstract. We discuss a new simple field theory approach of Coulomb systems. 
Using a description in terms of fields, we introduce in a new way the statistical 
degrees of freedom in relation with the quantum mechanics. We show on a 
series of examples that these fundamental entropic effects can help account for 
physical phenomena in relation with Coulomb systems whether symmetric or also 
asymmetric in valence. On the overall, this gives a new understanding of these 
systems. 



1. Introduction 

Recently, we have developped a simple field theory (FT) approach to describe liquids 
[1] and more specifically electrolytes at interfaces. The use of density fields as basic 
elements of description rather than particle density distributions comes together with 
a radically different way of accounting for the statistical degrees of freedom of the 
system. It is the quantum mechanics which sets the correct combinatory of the states 
by telling us which two states can be distinguished in the phase space. In this paper, 
on a series of examples related to Coulomb systems, we illustrate how within our 
approach the denombration of the states can help us interpret and give an original 
insight into a series phenomena related to Coulomb systems. In the first part of the 
paper, we discuss general aspects and consequences of the formalism. In the second 
part, we illustrate the entropic role of the quantum mechanics in Coulomb systems. 

2. A new and simple field theory to describe liquids, exact results 

The fundamental choice of our FT is to take directly the density of matter p(r) 
as the basic real valued field variable. As the densities fluctuate, we then fix 
the chemical potential /i and set the partition function to be the grand potential. 
It is then a functional integral according to —[3PV — ln2 = ln[/ Vp e~^^^''^]^ 
where j3 is the inverse temperature. The Hamiltonian functional H corresponds 
to a simple meaningful mean field approximation for the chemical potential /3/i = 
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ln(p(r)A3) + /p(r')w(|r - r'|)dr' which leads to l3H[p] = J p{r) [ln(p(r)/p^) - 1] dr + 
/ p{r)p{r')v{\r — r'\)drdr' where Pz = e^^^ /A? is the activity. In [T], we have shown 
that provided a renormahzation, this theory reproduces exactly the standard statistical 
mechanics. The peculiarity of the approach lies in the first term which is entropic and 
accounts for the quantum mechanics denombration of states with i) the elementary 
volume in the phase space A'' and ii) the principle of indiscernibility for particles 
both properties which are, for N particles, traditionally given by the J|-^ dr.i/ {N\K^) 
measure. 

Dyson like relations or equations of movement. A direct consequence of the 
integral functional formalism, can be seen in the invariance of the functional integral 
with respect to a change of the dummy field variable, for any functional A[p\ we have 
for the averages ((5 A = {A[p\{5[3H[p\/ 5p)) . Starting from the partition function, 
the application of this relation gives a new expression of the chemical potential: 
(3p, = (ln[p(r)A'^]) + (3 J v{\r — r'|)(p(r'))c?r'. This relation has been discussed in view 
of more standard expressions [2] . Where field theory shows a simple mean interaction 
potential term and where the average of the entropic logarithmic term includes more 
intricate correlations, in contrast, standard expressions emphasize the interaction 
potential, the logarithmic term being the trivial ideal system. This displacement 
in the way correlations are treated, naturally suggests new approximations. 

New perspectives on old exact relations. We have rederived in |2j the exact 
density contact theorem and the Virial theorem. We have pointed out that our 
approach provides yet another point of view on this relation emphasizing the role 
of the logarithmic functional in the demonstration. 

3. Entropic effects in Coulomb systems. 

Coulomb systems are peculiar mixtures: by this we mean that we have two species 
bound by the condition of electroneutrality. On the one hand, we have entropic 
properties related with densities of each of the two indiscernable species. On the 
other hand, in Coulomb systems we like to focus on the difference of the densities, the 
charge. Our point is to treat on the same footing the purely entropic properties and 
the one related to the potential as they appear at the same level in our Hamiltonian. 
And to show that this will lead to strong correlations between the different properties. 
Intuitively, when wc change the charge, we also modify the relative number of the 
species therefore their purely entropic combinatorics contribution. We hereafter 
discuss simple point ions to focus on these entropic terms. 

At a charged interface the mean field approximation gives the Gouy-Chapman 
theory. Discussing fluctuations beyond the mean field theory correspond to an 
expansion of the entropic part of the Hamiltonian which for Coulomb systems is 
jjent = Jdr {p+(r) [ln(p+(r)/p+) - 1] + p_(r) [ln(/9_ (r)/p_) - 1]} with separate 
terms corresponding to each ion. But for the more natural variables for coulombic 
systems, the charge q = p+ — p- and the total density s = + p^, the entropic 
functional mixes the fields as we have p± ~ (s ± 9)72. Note a peculiarity in 
FT is that linear combinations of the fluctuating fields q, s are natural and not 
simply combinations for the average quantities. In the expansion, the quadratic 
part of the Hamiltonian includes the Coulomb potential interaction. With this 
quadratic term, we obtain the Debye limiting law in the bulk. The specificity of 
the FT lies in the higher order terms which include products of q and s fields 
gffent _ —i/(^2p) J (7^(r)s(r) dr + ... in local couphng terms. 
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3.1. Application to inhomogeneous systems. 

Ionic depletion at a neutral interface^ For a neutral interface and a symmetric 
electrolyte, there is intuitively no profile for the charge at a neutral interface. However, 
the loop expansion of the theory shows a profile of the quadratic fluctuations of the 
charge with a simple expression with the distance to the wall and a dependance on 
rj — K'q/{87tp) which characterizes the Coulomb interaction strength where Kd is 
the inverse Debye length. Physically, some fluctuations are frustrated at the vicinity 
of the interface due to the absence of ions in the other half space. The role of the 
entropic term is to imply that this coulombic effect is coupled to the total density field 
although this quantity is not directly related to Coulomb interaction. The pertinence 
of this profile is further justified as it verifies the contact theorem which is not the case 
of the intuitive Gouy-Chapman or MSA (Mean Spherical Approximation) approaches. 

The anomalous capacitance behaviour. Experimentally and in numerical 
simulations, we find systems where the electric capacitance decreases at low reduced 
temperature. This effect is non intuitive as one can expect the electric response of 
the system to decrease with increasing thermal agitation which is what the Gouy- 
Chapman and MSA theories predict. This phenomenon appears at low reduced 
temperature when the Coulomb interaction becomes stronger. Within the field theory, 
the phenomenon can easily be understood with the previously described depletion 
in the density profiles which naturally accounts for the decrease in the electric 
capacitance. For Point ions [4], we have obtained a simple analytic expression of 
the capacitance. This has been corrected to account for the size of the ions [5]. The 
comparison with the numerical simulations is rather good and the simplicity of the 
expressions and of the corrections indicate we capture the significant physical effect. 

3. 2. More entropic effects for multivalency 

Entropic effects are also important for valence asymmetric systems, z+ : z_ electrolytes 
where the electroneutrality condition has now modified the number balance of the ions. 

Anomalous capacitance for valence asymmetry. In [6], we have considered 
capacitance curves for different valencies, the curves are rather scattered but one 
finds that one should first redefine the temperature T* T* = T* / {z^zJ) which 
is also a scaling with the ionic strength. However, the ionic strength z+z_ scaling 
does not discriminate between a 2:2 and a 4:1 electrolyte. The FT introduces a 
distinct parameter Zas = {z+ — z_ ) / ^/z+z- which is really characteristic of the 
valence asymmetry. As in simulation results, FT predicts an increase of the differential 
capacitance peak with Zag. This parameter appears naturally in FT as we can perform 
simple operations on the fluctuating fields and again it is related to the entropic term. 

Potential of zero charge (PZC) for size and valence asymmetric systems. For the 
neutral interface, numerical simulations results show the existence of a spontaneous 
polarisation due to the asymmetry in size and/or in valence The first is rather 
intuitive as the smaller ions by reaching nearer the interface induce a polarisation. Less 
intuitive is the polarisation due to valence asymmetry. Indeed, in figure 1, standard 
approaches in liquids state physics [7 do not account for such effect. On the contrary, 
FT predicts the polarisation due to multivalency and we show results obtained by 
adding the MSA theory for size asymmetry to the FT result for point ions (size 
asymmetry being immaterial for point ions) [8J. The sign of the PZC is in agreement 
with the charge contact theorem 'Q] and with the simulation results. The simplicity 
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Figure 1. PZC for 3:1 electrolytes as a function of d\/d2 ionic diameters ratio. 
The symbols are for the numerical simulation with the corresponding FT plus 
MSA predictions (see text), concentrations are 0.033, 0.333, 0.666 Mol/1 from top 
down. The dashed line are for the MSA in the reverse order for the concentrations. 

of the correction suggest again that we focus on the adequate physical effect. 

Ionic criticality and valence asymmetry. We have derived in the FT an equation 
of state for ionic systems. In the loop expansion, the pressure is a series in terms of ly. 
(3P = p (l — iyy — d^z^g-q^ + ...). The linear term —77/3 is the standard Debye term, 
further we have a simple dependence with the Zas parameter. We have generalized 
this expression to account for hard sphere effects by using the MSA pressured and 
calculated the asymmetry contribution do using the MSA correlation function. The 
term in the pressure associated with valence asymmetric systems gives a simple 
description of the decreasing critical temperature and increasing critical density in 
agreement with Monte Carlo (MC) results [TO]. Respectively for the 1:1, 2:1, and 
3:1 electrolytes the critical temperatures and densities are T* — 0.079, 0.060, 0.038 
and p* = 0.0145,0.021,0.034. This seems to indicate that the physics is likely to be 
associated with the simple number balance effect found in multivalent systems. 

4. Conclusion 

Within a FT framework, we have shown how the formalism can introduce a new point 
of view on Coulomb systems. The basic and fundamental quantum mechanical count- 
ing of states can lead to important consequences in Coulomb systems which have at 
least two indiscernable species. These entropic effects lead to a unavoidable correlation 
between charge and density fluctuations. This simple statement is sufficient to give 
meaningful interpretations, we believe this is useful to complete the standard liquid 
state approaches generally focused on excluded volume effects. 
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§ Corresponds to replacing Kjy — > 2r the MSA screening distance parameter, 2r = ^1 + 2Kj:) — 1 



